Characterizing microstructure of living tissues with time-dependent diffusion 
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Molecular diffusion measurements are widely used to probe microstructure in materials and living organisms 
noninvasively The precise relation of diffusion metrics to microstructure remains a major challenge: In complex 
samples, it is often unclear which structural features are most relevant and can be quantified. Here we classify the 
structural complexity in terms of the long time tail exponent in the molecular velocity autocorrelation function. 
The specific values of the dynamical exponent let us identify the relevant tissue microanatomy affecting water 
diffusion measured with MRI in muscles and in brain, and the microstructural changes in ischemic stroke. Our 
framework presents a systematic way to identify the most relevant part of structural complexity using transport 
measured with a variety of techniques. 



When modeling or interpreting bulk transport in realistic 
disordered samples, e.g. composites, porous rocks, or living 
tissues, the challenge is to identify what part of microstructure 
to focus on [1-8]. Macroscopic transport is affected by multi- 
ple aspects of the immense microscopic complexity, yet their 
relative importance is hard to estimate and compare. Here, we 
describe how a measurement itself may give us a hint regard- 
ing which parts of the microstructure are most relevant and, 
thereby, can be quantified. To that end, we suggest to employ 
the long temporal correlations in molecular diffusion which 
preserve the footprint of the underlying structural complexity. 
These correlations manifest themselves in the power law tail 
of the molecular velocity autocorrelation function 



V{t) = {v{t)v{0)) ^ t-^^^^\ ^>0. 



(1) 



Practically, the tail (1) can be identified in the way the time- 
dependent instantaneous diffusion coefficient 



Anst(t) 



d {5x^) 
dt 2 



Jo 



) Doo^const-t-"^ (2) 



approaches the finite bulk diffusion constant I^oo • The quan- 
tity I^inst(^) is accessible with any technique [3, 8] measuring 
the mean square molecular displacement {5x^{t)) in a partic- 
ular direction, see Eqs. (4)-(6) in the Supplemental Material 
[9]. 

Structural complexity (disorder) presents itself in many dif- 
ferent forms, e.g. Figs. 1 and 2. Our key result is the relation 



^ = (p + d)/2 



(3) 



(see [9]) between the dynamical exponent ^ in equations (1) 
and (2), and the structural exponent p characterizing global 
structural organization in d spatial dimensions. 

The structural exponent p determines the T{k)\k^Q ~ k'^ 
behavior of the Fourier transform of the correlation function 
r(r) for the underlying microstructure. Hence, p character- 
izes global structural complexity, taking discrete values ro- 
bust to local perturbations. This enables the classification of 
the types and topologies of the disorder. Figs. 1 and 2. 



The relation (3) provides a way to determine the exponent 
p and, thereby, the structural complexity class, using any type 
of bulk diffusion measurement. Local properties affect the co- 
efficients, e.g. the values of I^oo and of the pref actor of 
in (2), but not the exponent i?. The latter is robust with respect 
to variations between samples of a similar origin. This pic- 
ture is akin to critical phenomena [10], where the phase tran- 
sition temperature is non-universal (sensitive to microscopic 
details), whereas the critical exponents distinguish, based on 
global symmetries, between universality classes of long-range 
fluctuations. Here, we describe how transport distinguishes 
between the classes of stationary structural fluctuations. 

Figure 1 illustrates how diffusion distinguishes between the 
disorder classes via the relation (3) in d = 1 dimension. The 
Monte Carlo simulated diffusion is hindered by the perme- 
able barriers with mean density n and permeability n. The 
structural complexity classes, embodied in the types of bar- 
rier placement (Fig. la), exhibit qualitatively different struc- 
tural exponent p in the barrier density correlator (Fig. lb). 

Order (periodic placement in Fig. la) is reflected in the 
Bragg peaks in r(/c), with F = for /c below the minimal re- 
ciprocal lattice vector, formally corresponding to p = oo. As 
coarse-graining beyond the largest lattice vector does not con- 
tribute to the structural fluctuations, V(t) decays and I^inst(^) 
reaches I^oo exponentially fast (^ = oo) already at the lattice 
scale, with the decay rate determined by the lattice specifics. 

Structural disorder can be introduced in qualitatively differ- 
ent ways. Hyperuniform disorder [11, 12] is characterized by 
suppressed long-range fluctuations, with the variance in the 
number of restrictions in a given domain increasing slower 
than the domain volume (sub-Poissonian statistics), reflected 
in p > 0, and ^ > d/2. Here, we displaced the barriers 
independently from their positions in a regular lattice by in- 
dependent random shifts leading to V{k) ~ /c^ for small k, 
such thatp = 2 (Fig. lb), yielding = 3/2 (Fig. Ic) accord- 
ing to equation (3). Short-range disorder corresponds to the 
variance of the number of restrictions scaling with the mean 
number in a given domain, consistent with the central limit 
theorem. Here, we chose each successive interval am be- 
tween barriers independently from the distribution P{a) with 
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FIG. 1. Time-dependent diffusion distinguishes between structural 
complexity classes in one dimension, represented by the placement 
of identical permeable barriers with the same mean density. (A) 
Order (red), hyperuniform disorder (green), short-range disorder 
(blue), and strong disorder (magenta). (B) The barrier densities have 
qualitatively different large-scale fluctuations, reflected in the small- 
k behavior of their density correlator r(A:) ^ (see text). (C) 
Numerical results confirming the relation (3). The time-dependence 
(2) clearly distinguishes between the four arrangements, while the 
value i^oo is the same for all of them. The dashed lines are the exact 
power laws from equation (9), and the exponential decrease is from 
the exact solution (see text and [9]); Tr = a /2k. 



mean a = l/fi and finite variance a^. This results in the fi- 
nite plateau r|/e^o = cr^/a^, such that p = 0, qualitatively 
similar to Poissonian disorder. Hence, ^ = 1/2 (Fig. Ic), 
and, generally, ^ = d/2'm d dimensions [13, 14], cf. equa- 
tion (3). Finally, strong disorder, with structural fluctuations 
growing faster with volume than prescribed by the central 
limit theorem, is reflected in a diverging r|/c^o, i-e. the ex- 
ponent p < 0, and li) < d/2 (weak self-averaging). Here, 
we used the Levy (fat tail) distribution P{a) ~ with 
^ = 7/4 for the successive barrier intervals, such that the vari- 
ance ((a — a)^)p diverges. This yields p = /i — 2 = — 1/4 
and ^ = (/i — l)/2 = 3/8 in agreement with equation (3). 

Higher dimensions d > 1 provide more ways to realize 




FIG. 2. Extended structural disorder classes, d > 1. (A) Randomly 
placed and oriented permeable barriers (lines), = 1 in d = 2 
dimensions. (B) Randomly placed and oriented rods, = 1 in d = 
3. (C) Randomly placed and oriented permeable barriers (planes), 
^5 = 2 in d = 3. (D) Structure correlator r(A:) ^ (numerically 
calculated and angular averaged, arb. units) for (A) (magenta), (B) 
(grey), and (C) (green) exhibits the negative exponent p = —ds. 



the same basic disorder classes. Various periodic arrange- 
ments would yield the same qualitative behavior, p = oo and 
"3 = oo. Hyperuniform disorder can be realized for different 
p > 0. While p = 2 when the restrictions are independently 
displaced away from the lattice sites, p c:^ 1 for a maximally 
random jammed state [12]. Equation (3) provides the possi- 
bility to observe the jamming transition, fromp = to p = 1, 
via diffusion in-between packed impermeable beads. 

Remarkably, higher dimensions open up ways to realize 
strong structural fluctuations, with diverging T{k), corre- 
sponding to p < and i9 < d/2, without a need to invoke a 
Levy distribution. A negative p (Fig. 2) can be achieved very 
naturally, by organizing micro structure in terms of randomly 
placed and oriented regular components (e.g. infinite lines, 
planes) with dimensionality ds < d, in which casep = — (a 
negative integer), and 2i? corresponding to their co-dimension. 
The first such example [15] is the extended disorder realized 
by random permeable hyperplanes, ds = d — 1, Fig. 2a, re- 
sulting in ^ = 1 /2 in any d. Likewise, randomly placed and 
oriented rods, ds = d — 2, embedded in d = 3 dimensions 
(Fig. 2c) would realize p = —1 and = 1 < 3/2. The above 
examples merely represent each disorder class; e.g. the "rods" 
from Fig. 2c can be structurally complex, permeable or imper- 
meable. What matters is the long-range correlations. Clearly, 
for components with finite extent, the disorder becomes short- 
ranged, i9 d/2, when the rms molecular displacement ex- 
ceeds their size. 

Above, we assumed that the molecules (the random walk- 
ers) can spread everywhere. When impermeable boundaries 
split the space into disconnected parts, equation (3) applies 
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FIG. 3. Time-dependent diffusion transverse to muscle fibers [16] 
(d = 2) reveals extended structural disorder class of = 1, pro- 
vided by the muscle fiber membrane (sarcolemma). (A) The longi- 
tudinal, Ai, and the transverse, Ax, diffusion tensor components for 
calf tongue genioglossus (TG, blue circles) and heart (H, red dia- 
monds). Solid lines are the fit of X±{t) to Supplementary Eq. (30) 
with d = 2. For fit results see Table 1 [9]. (B) Data for A±(t) re- 
plotted as function of t~^/^ consistent with i9 = 1/2. Equation (3) 
yields p = — 1; hence, ds = 1 (see text and Fig. 2d). (C) Histo- 
logical slice [17] of skeletal muscle across the fibers. Note the tight 
cell packing achieved by straight cell walls, qualitatively similar to 
Fig. 2a. 



separately to the contribution from each part, which then add 
up. The most relevant disorder contribution is the one with the 
smallest i?, yielding the slowest power law tails (1) and (2). 

As a result, measuring the exponent ^ with any time- 
dependent diffusion technique allows one to determine the 
disorder universality class via the exponent p using the rela- 
tion (3). Let us now apply this framework to diffusion mea- 
sured with IMRI [3] (dlMRI) in tissues. The dlMRI in muscles, 
Fig. 3, reveals strong disorder (weak self-averaging) in d = 2 
with ds = 1, realized by the sarcolemma; dlMRI in the brain. 
Fig. 4, reveals short-range disorder along the d = 1 neurites. 

In Fig. 3, we analyze the time-dependence of diffusion ten- 
sor eigenvalues in the fresh ex vivo muscle tissue samples 
measured by Kim et al [16]. The nondispersive eigenvalues 
Ai correspond to the unrestricted diffusion along the fibers. 
The transverse components X±{t) in the two-dimensional 
fiber cross-section. Fig. 3c, are strongly dispersive. Repre- 
senting the data as function of we observe the asymp- 
totic tail (2). Indeed, the fit of X_\_ {t) to Supplementary Eq. (5) 
yields ^ 0.5 for both tongue and heart (Fig. 6 [9]), exempli- 
fying weak self-averaging, i9 < d/2, in contrast to = 1 ex- 
pected for the d = 2 short-range disorder. We thus conclude 
that the restrictions to water diffusion are strongly spatially 



FIG. 4. Dispersive diffusion in cerebral gray matter reveals strong 
structural correlations. Notably, they remain qualitatively the same 
after global ischemia. (A) Original data [20] for 'd-sin' and 'cos' 
gradient waveforms, fitted to Supplementary Eq. (6), yields i9 = 
0.50 ± 0.07 for normal and i9 = 0.49 =b 0.05 for postmortem 
brain, consistent with the p = short-range disorder along the one- 
dimensional neurites (dendrites and axons) in both cases. The role 
of the disorder (the slope) increases after ischemia. (B) Varicose ax- 
ons from rat hippocampus area CAl, ref. 22, rationalizing the d = 1 
diffusion inside narrow randomly oriented disordered neurites. 



correlated on the scale of the diffusion length (up to ~ 30 /im 
in this measurement), which puts them into the extended dis- 
order class of Fig. 2a with ds = I. 

In [9], we argue that the relevant restrictions are in fact 
muscle cell membranes (sarcolemma), and quantify their per- 
meability and cell size (Table 1 [9]). The good agreement 
between the fit parameters and histological values can be ra- 
tionalized by comparing a typical histological slice transverse 
to muscle fibers [17] (Fig. 3c) with the random barriers in two 
dimensions (Fig. 2a). Tight packing of muscle cells makes the 
fiber walls fairly flat and spatially correlated even over length 
scales exceeding typical fiber diameter, qualitatively affect- 
ing the structural correlations, r{k) ~ 1/k, within the plane 
transverse to the fibers. 

Hence, the dynamical exponent (1) establishes the effect of 
cell walls on the IMRI signal, leading to the first non-invasive 
method to image cell membrane permeability in vivo, with the 
potential to correlate with tissue physiology and pathology. 

We now turn our focus to brain. Fig. 4. The early observa- 
tion [18] of the almost two-fold decrease in the water diffusion 
coefficient minutes after brain injury has helped to spur the 
development of in vivo dlMRI. IMeasured at times t ~ 100 ms 
when any residual time dependence in brain is small, the dif- 
fusion coefficient is now widely utilized clinically as a non- 
invasive diagnostic marker for acute ischemia [19]. However, 
the microscopic origin of this phenomenon has remained un- 
der debate for two decades. A closely related challenge is 
to identify the predominant restrictions or cell mechanisms 
which determine water diffusion in healthy brain. 

Here, we address these two related questions by focussing 
on the Fourier transform V{uj) of the velocity autocorrelation 
function (1) [Supplementary (6)] in rat cortical gray matter. 
We observe that the real part of V{uj) measured by Does et al. 
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[20] with oscillating gradients [3] exhibits the i9 = 1/2 disper- 
sion (Fig. 4) in the whole frequency range, cj/27r < 0.5 kHz. 
This value of "3 is striking for two reasons. First, microstruc- 
ture must be highly correlated, since, for the isotropic [20] dif- 
fusion ind = 3, the naively expected ^ = 3/2 in the absence 
of structural correlations, p = 0. Second, the value "3=1/2 
is the same before and after global ischemia. 

In [9], we radically narrow down the scope of plausible sce- 
narios for the observed change in V{uj) based on the previ- 
ously unidentified "0 = 1/2 tail (1), ruling out active stream- 
ing breakdown, full confinement by impermeable walls, cell 
swelling, and increase in the cytoplasmic viscosity as the pre- 
dominant mechanisms. Out of the remaining possibilities, we 
argue that most of the dispersion arises from the diffusion 
along randomly oriented narrow neurites (mostly dendrites 
in gray matter, and, possibly, some axons). Strong apparent 
correlations are maintained by the impermeable neurite walls 
such that, effectively, d = 1, while the disorder along the neu- 
rites (such as shown in Fig. 4b) is short-ranged, p = 0, so that 
the universality class is that of Fig. 1 (blue). 

This disorder, for the dendrites, may include [21] spines, 
variations in thickness ("beads"), and in local directionality 
on the ^ 1 /im scale; for the axons, the synaptic boutons (vari- 
cosities) separated [22] by 3 — 6 /im. Our short-range disorder 
conclusion is remarkably consistent with the measured vari- 
ance in the varicosity number within a window growing in 
proportion to the mean within this window, ref . 22, a defining 
signature of the p = exponent. Ischemia causes beading, 
i.e. more pronounced varicosities in both dendrites and ax- 
ons [23, 24], which is likely to increase the disorder. This is 
consistent with the increase in the pref actor in the cj^/^ con- 
tribution to T>{uj), Fig. 4, as this pref actor generally grows 
when disorder gets stronger. Our analysis yields that the effec- 
tive "permeability" of shafts between beads drops more than 
twofold in ischemia. 

Here, our approach underscores the value of the time de- 
pendence of diffusion, rather than of a single number = 
V{uj)\^^=o, for uncovering the origin of a complex biophysi- 
cal phenomenon. It adds a crucial piece, the short-range dis- 
order, to the picture of impermeable cylinders for the neurites 
[25, 26], and is consistent with the decrease in Doo and ap- 
pearance of beads under a mechanical stress in ex vivo axons 
[27]. The present framework may stimulate more focussed 
investigations of ischemic stroke, as well as of other neuro- 
logical disorders. In particular, one could correlate the time- 
dependent diffusion with the morphological changes during 
status epilepticus and electrical activation [28], and severe hy- 
poglycemia [29], also known to reduce the value of Doo- The 
reduction in the number of axonal varicosities in Alzheimer's 
disease relative to the healthy brain [30] is likely to result in 
the decrease of the o;^/^ contribution due to the reduced struc- 
tural disorder. 

To conclude, we connected the dynamical exponent (3) 
to the global structural organization, in order to study mi- 
crostructure with any type of diffusion measurement. This 
framework is particularly useful for biological tissues: While 



biophysical parameters may vary strongly and continuously 
between samples, the exponent ^ takes fixed values deter- 
mined by the disorder universality class, and is robust with 
respect to the biological variability. As a result, we iden- 
tified the dominant role of cell membranes restricting water 
motion in muscles, and argued for an increase in the structural 
disorder along the neurites as a cause of the diffusion coeffi- 
cient decrease after ischemic stroke. We believe the presented 
classification of the disorder could help identify and quantify 
the dominant types of restrictions in other living tissues, as 
well as in classical diffusion or heat or electrical conduction 
in composite materials, porous media, and other structurally 
complex samples. Extending this approach to the quantum or 
wave transport would tie the rich physics of localization with 
the types of the global structural organization. 
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SUPPLEMENTARY MATERIAL 

Time-dependent diffusion from structural disorder 

Diffusion metrics. The fundamental quantity, the velocity 
autocorrelator (1), is often difficult to measure directly. In- 
stead, there exist a number of equivalent time- or frequency- 
dependent diffusion metrics, with the relations between them 
described in detail in refs. 31 and 32. To interpret vari- 
ous kinds of diffusion measurements, such as dMRI results 
[16, 20], here we outline how the power law tail (1) manifests 
itself in these metrics. From the outset, we assume the sam- 
ple to be statistically isotropic, so that the diffusion metrics 
are isotropic tensors, and the correlation functions depend on 
r = |r| and k = |k|. Generalization to the anisotropic case 
presents no conceptual difficulty, but makes the presentation 
more cumbersome. 

The instantaneous diffusion coefficient I^inst(^) defined in 
equation (2) of the main text is the natural metric to study 
structural correlations, as it quantifies how the spreading of a 
packet of random walkers is hindered by the microstructure at 
the time scale t. From our perspective, it is a perfect quantity 
to determine the exponent ^. However, this is not the most 
commonly utilized metric in practice. 

The most commonly reported diffusion coefficient 

D{t)^&^ = jJ^D,Ut')dt' (4) 

describes the dynamics of the cumulative, rather than instan- 
taneous, mean squared displacement along a particular direc- 
tion X over the diffusion time t. This is the case both in the 
dMRI [3, 16] and in the direct molecular tracking techniques 
[8]. This definition has a perceived advantage of dividing by 
time, rather than differentiating with respect to it: clearly, dif- 
ferentiating increases the noise, while dividing does not. 

However, the definition (4) may mask the exponent i^. In- 
deed, defined in this way, the long t behavior 

D{t) -D^^ const • t"^ , ^ = min 1} . (5) 

In other words, for the tail (1) to be manifest in D{t), it should 
be sufficiently slow, < 1; in this case it "goes through" the 
averaging over the increasing interval t in equation (4). In 
the opposite case, ^ > 1, the term in I^inst(^) becomes 
subleading to the 1/t term from the integral in equation (4) 
converging at short t. 

Hence, to practically determine the dynamical exponent 
one could first check whether the fit to equation (5), using 
the "less noisy" definition (4), produces the value i9 < 1. If 
it does (as in our example of diffusion transverse to muscle 
fibers), this is it, ^ = i^. In the opposite case, the fit would 
yield the 1/t tail, ^ = 1, which would mask the true value of 
^ > 1. Then, one must perform the differentiation I^inst(^) = 
dt[tD{t)] and obtain from the fit to equation (2), with the 
unfortunate effect of amplifying the measurement noise, as 
shown by comparing Figs. 1 and 5. Practically, this results in 




FIG. 5. Cumulative diffusion coefficient, equation (4), for the one- 
dimensional example of Fig. 1. Dashed lines correspond to the 
asymptotic power law decrease of D{t). For ^ = 1/2 and = 3/8 
(short-range disorder, blue, and strong disorder, magenta), the power 
law in D{t) coincides with that in Dinst (t), cf. equation (5), whereas 
fori9> 1 (periodic, red, and hyperuniform, green), it is masked by 
the 1/t term. Taking the derivative Dinst = dt{tD(t)) reveals the 
values of ^ (cf. Fig. 1) but increases the noise. 

more stringent requirements on the signal-to-noise ratio and 
on the greater number of experimental time points. 

There is another useful way of uncovering the exponent ^, 
as long as 1? < 2, without the need to take a time derivative. 
The same power law tail 

RtV{uj) = \{v-^v^) ^ Doo^con^i'\ujf , ^ , (6) 

persists in the dispersive diffusivity V{(jj) = j^^t e*^^P(t), 
which is the Fourier transform of the retarded velocity auto- 
correlator (1). The physical meaning of V{(jj) is in relating the 
current Ju;,r = — I^(cj) VrV^u;;r of the random walkers to their 
density gradient, somewhat similar to the dispersive electrical 
conductivity; it defines the pole of the diffusion propagator, 
see refs. 13-15, 31, and 32 and also the discussion below. Re- 
markably, there exists a standard dMRI measurement proto- 
col, the oscillating gradient technique [3, 20], which directly 
measures [32] ReP(a;). This is the quantity used in the ex- 
ample of diffusion in cerebral gray matter. Fig. 4. 

Derivation of equation (3). Homogenization. In this 
work, we consider the most widespread situation, when a 
sample has a nonzero macroscopic diffusion constant I^oo = 
[{5x^) /2t]t^oo^ i-e. the diffusion asymptotically becomes 
normal, or Gaussian. A well-defined macroscopic I^oo, ob- 
served in an overwhelmingly broad variety of microscopically 
heterogeneous samples, attests to the robustness of the dif- 
fusion as a Gaussian fixed point with respect to adding the 
structural complexity (disorder). In this case, a macroscopic 
sample represents the disorder ensemble, i.e. the system is 
self-averaging [33]. Conversely, the absence of I^oo, e.g. for 
fractals, near a percolation threshold [2, 34], or for random 
drifts in one dimension [35, 36], signifies so-called anoma- 
lous diffusion [2] . 

The general relation of the long-time behavior (1) and 



2 



(2) to the micro structure rests on the homogenization ar- 
gument: At long diffusion time t, the sample, as seen by 
random walkers traveling over a growing diffusion length 
L{t) = {5x^{t))^^'^ 0^ \^2Doot, appears increasingly more 
uniform due to self- averaging. The sample is being effec- 
tively coarse-grained over L{t), such that the strong micro- 
scopic heterogeneity is gradually forgotten, and the deviation 
^i^(r) = i^(r) — Doo of the smoothly varying coarse-grained 
diffusion coefficient D{r) from Doo becomes small. This jus- 
tifies calculating the self-energy part of the disorder-averaged 
diffusion propagator only to the lowest (second) order [31] in 
the variable component SD{r). Eventually, the perturbative 
treatment around D^o becomes asymptotically exact, and the 
residual deviation 
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TD{k) 
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is given in terms of the Fourier transform Voik) = 
jd^Y e~^^^V d{t) of the two-point correlation function 
^D{r) = {5D{vo-\-y)5D{yo)) in spatial dimensions. Using 
the relation between V{uj) and D{^^i{t), 
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(which can be derived using the cumulant expansion as out- 
lined in ref. 31), we obtain 



AnstW-i^o 



dDoo J 



d^k 



TD(k)e 
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Equivalently, the latter can be recast in the form 



Anst (t) ^ i^oo + const • {{5Df)\L{t) , (10) 
where {{5DY)\l is the variance of the Gaussian-smoothed 

values (5I)(r)|L = /dV (5I)(r + r') e-^''/^V(^^')^^'- In 
other words, the diffusion effectively applies a low-pass fil- 
ter ^ 1^ to the Fourier components of D{y) and, thus, 
to its correlator Vjy(}i), admitting harmonics with progres- 
sively smaller wavenumbers k < 1/L{t). As the variance 
{{SD)'^)\l ^ decreases due to the smoothing, the mea- 
sured diffusion coefficient Dinst (t) monotonically decreases 
towards Doo- The power law exponent (3) is then directly 
related to the dimensionality d and the exponent p which de- 
termines the A: ^ behavior of T oik) ~ k^. 

We are interested in the spatial correlations r(r) = (n(ro + 
r)n(ro)) of the underlying micro structure n(r) responsible 
for the heterogeneity of D{r). Depending on the sample, n(r) 
may stand for the density of grains, traps, barriers and other 
structural components (e.g. Figs. 1 and 2), which is often 
strongly heterogeneous at the microscopic scale. Certainly, 
the coarse-grained D{r) is not equal to the local average of 
the strongly varying microscopic diffusion coefficient caused 
by n(r). However, the statistics of the large scale fluctuations 
of n(r) asymptotically approaches that of the coarse-grained 
D{r), such that for k ^ 0, 

TDik) C{n) . T{k) , C{n) = {dDp^/dnf . (11) 



This asymptotically local relation rests on the self- averaging 
property which ensures the smooth dependence Doo (^) on the 
sample mean n = (n(r)) of the restrictions. Hence, after 
coarse-graining, a typical small local fluctuation 5D{y) be- 
comes asymptotically proportional to the typical small local 
fluctuation n(r) — n, as long as the self-averaging assump- 
tion holds. (Conversely, singular dependence Doo{ti), e.g. at 
the percolation threshold, is associated with the lack of self- 
averaging.) In this way, the exponent p characterizes long 
range correlations in sample's micro structure, T{k) ~ /c^, and 
becomes accessible with a time-dependent diffusion measure- 
ment via the relation (3). 

Strong self-averaging in d dimensions occurs when the vari- 
ance {[5DY)\l decreases as the inverse "diffusion volume" 
L~^, such as for the short-ranged disorder (jp = 0), or faster, 
as for order or hyperuniform disorder (p > 0). Weak self- 
averaging corresponds to the decrease ^ L"^^ with a smaller 
power, < 'd < d/2. For p < —d, very strong fluctua- 
tions destroy self- averaging, a large sample does not repre- 
sent a disorder ensemble, the macroscopic Doo is undefined, 
and the present approach fails. Diffusion becomes anoma- 
lous [2], with mean squared displacement {5x^{t)) ~ t^/^ as 
t ^ oo characterized by the dynamical exponent 2; 7^ 2, see 
e.g. refs. 34-36. 

One-dimensional disorder universality classes. Here we 
provide the details of the analytical and numerical calculations 
used to obtain the results in Figs. 1 and 5. 

Monte Carlo dynamics. — For each disorder class, Monte 
Carlo (MC) simulated random walks of 4 x 10^ random walk- 
ers evenly split between 40 disorder realizations of TV = 1000 
barriers each, were used to average {Sx'^) over the paths and 
over the ensemble. The total length of each disorder realiza- 
tion could be either smaller or larger than TV a since barrier 
intervals were random (as described below). The trajectory of 
each random walker was a sequence of moves in a randomly 
chosen direction over a distance dx = y/2Dodt = 0.008 a 
during a time step dt, where Z^o is the unrestricted (free) dif- 
fusion constant. This choice of dx <C a ensured that the free 
diffusion was well simulated within each inter-barrier interval, 
i.e. the effects of the finite step dx at the scale of inter-barrier 
separation were already negligible. In this way, the added bar- 
riers can be viewed as the restrictions, or the "disorder", for 
the ideal free one-dimensional diffusion. 

The barrier permeability k (the dimensions of velocity) de- 
termines the dimensionless disorder strength, [15] 



(12) 



The value C = 1 was chosen for all disorder realizations 
and types. Keeping the same a, Do and n for all disor- 
der classes yields the same macroscopic diffusion constant 
Doo = ^0/(1 + C) = Do/2 for all simulations. 

Transmission across a barrier occurred with probability 
P (X ndx/D^ <C 1 (ref. 37). The total diffusion time was 
lOOr-p, corresponding to a maximum of 1.5625 x 10^ time 
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steps per walker, where 



density rik, entering the correlator (15), 



V 



a 

2^ 



(13) 



is the mean residence time within an average inter-barrier in- 
terval characterized by its surface-to-volume ratio S/V = 
2/a. The MC results were carefully calibrated to yield the 
exact result for Doo (with better than 1% accuracy), as well as 
using the exact result for Dinst(0 the periodic barriers, cf. 
equation (29) below and Fig. Ic, red curve. 

The random walk simulator was developed in C-i~i-. Simu- 
lations were carried on the NYU General Cluster, using 120 
central processing unit cores simultaneously, within a total 
time of about lOh per each disorder class. 

In order to obtain i^inst(^) in Fig. 1, the time derivative in 
equation (2) was calculated using the Savitzky-Golay smooth- 
ing procedure written in Matlab, with the 6th order polyno- 
mial interpolation over a window increasing with t to sup- 
press the MC noise that becomes relatively more pronounced 
at longer diffusion times. 

Asymptotic behavior of i^inst(^)- — Adding the barriers 
corresponds to the microscopic density of restrictions (dis- 
order) n{x) = Ylm^i ^{^ ~ ^m), where Xm are the barrier 
positions. In the Fourier domain, this results in the density 
variation Sn = n — n 
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The asymptotic analytical expressions (dashed power law 
lines in Fig. 1) are based on the relation (9) and on the "lo- 
cality" relation (11) asymptotically valid for A: ^ in any d. 
In our case. 



r(A:) = — {Sn-k^Uk) 



(15) 



is the barrier density correlator, and V is the system volume 
(length). Using equations (9), (11), and Doo{n) = Dq/{1 -\- 
nDo I k) valid in (i = 1 for any barrier placement, we obtain 
the exact expression for the small relative deviation from Doo 
in one dimension: 
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We will now use the general relation (16) with r{k) for the 
different disorder classes to obtain the corresponding asymp- 
totic behavior, dashed straight lines in Fig. Ic. 

Short-range disorder — A simple way to realize short-range 
disorder is to place barriers sequentially with their successive 
intervals am = ^m+i — i-i-d. (independent identically 
distributed) random variables chosen from a given probability 
density function (PDF) P{a) with a finite mean and variance. 
We now relate the limit r{k)\k^o, entering equation (16), to 
the parameters of P(a). The second moment of the barrier 
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with Xm — Xm' = ^^l^rn' avcragcd over the disorder 

by splitting the double sum, in the limit TV = nV oo, into 
three terms, with m = m\ m > m' and m < m\ where 
s = m — m\ and using the geometric series formula: 
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where Pk = /da e~*^^ Pis^) is the characteristic function of 
the PDF P(a), and * stands for the complex conjugation. The 
last term cures the uncertainty at /c = 0, setting r(A:)|/e=o = 
0, and does not affect the behavior in question at small but 
finite k. Combining the fractions, we find the barrier density 
correlator 
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PkPk 



(17) 



As the long-time diffusivity behavior depends only on the 
k^^ limit of T{k) (different from the value T{k)\k=Q = 0), 
we can represent this limit in terms of the mean and variance 
of any P{a) via its cumulants, pk = e~*^^~^ ^ /2+... 
stituting into equation (17) and taking the limit ^ 0, we 
obtain the limit r{k)\k^o = cr^/a^, as quoted in the main 
text. Substituting this limit into equation (16), we obtain 
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(dashed line on top of the blue MC line in Fig. 1). The param- 
eter cr^/a^ = 1/4 for our choice of the lognormal P(a). The 
usage of the lognormal distribution here is not important; any 
P{a) with finite variance and the same ratio a /a would yield 
the same asymptotic dependence (18). 

Strong disorder — When the interval PDF has a "fat tail" 
so that the variance a diverges, the above result (18) does not 
apply. The asymptotic form of the Levy distribution tail [2] 



P(a) 
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1 < /i < 2. 



(19) 



This choice for /i ensures the finite mean d = J^daaP{a) 
but yields an infinite variance due to the fat tail. Practically, 
we generated the disorder according to the interval PDF (19) 
by choosing a = amin where ^ is the standard random 

variable uniformly distributed in [0, 1], and amin = a(/i — 
l)//i, such that a > amin and C = /ittmm- 

Plugging the characteristic function [2] for a > 0, pk = 
^-^ka-c'k^[l+itM^^,/2)] ij^^Q equation (17), we find the k 

divergency in the correlator 
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(20) 
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In the second equation above, r(l + fi) is Euler's Gamma 
function. Equation (20) shows that, by using the Levy PDF 
(19) for the i.i.d. intervals, it is possible to generate a medium 
with an infrared-divergent density correlator. Of course, this 
is just one particular way of doing so, but it suffices here, as 
we need to describe only one member of the corresponding 
disorder universality class. 

The singular barrier correlator (20), substituted into equa- 
tions (7) and (16), causes the power law dispersion with the 
exponent -i^ = {li — l)/2. 
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Remarkably, the dispersion with the exponent < < 1/2 
decreases qualitatively slower than that for the finite vari- 
ance. This observation can allow one to determine the power 
11= 1 + 2^ of the Levy- stable PDF of intervals from the time- 
dependent diffusion. The "finite- variance" exponent ^ = 1/2 
is reached for the borderhne case of /i = 2, separating infi- 
nite and finite values of cr. In this limit, setting C cr^/2, 
equations (21) and (22) correspond exactly to equation (18). 

Finally, we can express the pref actor of in equation (22) 
solely in terms of the single barrier properties and the expo- 
nent 'd, by excluding the Levy tail normalization C\ 
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(23) 

Equation (23) for 7^ = 3/8 agrees well with the MC simula- 
tions in Fig. Ic. The cumulative D{t) — Doo, obtained from 
this equation using the relation (4) [which amounts to divid- 
ing by 1 — 7^], agrees well with the corresponding cumulative 
MC-generated diffusion coefficient in Fig. 5. 

Hyperuniform disorder. — We realize this disorder class in 
d = 1 by displacing the barriers from their positions in a pe- 
riodic arrangement. Now it is the displacements that are i.i.d. 
random variables taken from the PDF Pdispi(0- Clearly, the 
long-range order is preserved, as the system on average "re- 
members" about its ideal lattice positions. Hence, one expects 
the structural fluctuations to be qualitatively less pronounced 
as compared to the disorder types described above. 

To calculate the correlator (15) for the hyperuniform place- 
ment, we begin from the density (14), where now Xm = 
are the i.i.d. displacements generated ac- 
cording to Pdispi(0- Using 8{k)\k^Q = V/27T where V is the 
system length and I = N = nV, we obtain 

N 

m,m' = l 

^{27rnf^S{k) - 2 • 27rnNS{k) . (24) 

Averaging of the double sum over the disorder is done by split- 
ting it into the part with m = m' yielding N, and the part with 



m ^ m' yielding the Debye- Waller factor (e 



|Pdispi,/cP, where pdispi,/c = /d^e '^^Pdispi(0' multiplied 
by 
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The sums in the right-hand side span over all the reciprocal 
lattice vectors km = 27rm/a, m = 0,±1,±2,..., which is a 
consequence of the Poisson summation formula e*"^^^ = 
^ ^{k - km) valid in the limit N ^ 00. Putting all the 
pieces together, we obtain 
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m/0 



(25) 

This correlator is familiar from the X-ray scattering in crys- 
tals: a series of spikes with decreasing amplitude, together 
with the "incoherent" background determined by the dis- 
placement PDF. Using the cumulant form of Pdispi,/c = 
g-c^dispi^ we obtain the universal k ^ behavior 



Substituting (26) into equation (16), we finally obtain 



(26) 
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In our MC simulations, we took Pdispi = l/a for —a/2 < 
^ < a/2 and zero otherwise; substituting its variance cTdispi ~ 
a^/12 into equation (27), we obtain the dashed line agreeing 
with the green MC curve in Fig. 1. 

Order (periodic lattice). — The ordered (periodic) case ad- 
mits the exact solution [38, 39] for the diffusion propagator 
Guj,q in terms of the Bloch waves. In our present notation, 
with a = a the lattice period and ( = 2£/a = Dq/ kg, it reads 



G, 



uj,q 



t'o (1 + f ) + knacotkna {kl - q^ f + D^kl 

(28) 

where A;„ are the positive roots of the equation cos kjiCL — 
cosga + kn^^sinknCi. One readily checks that the n = 
term of the sum (28) yields the macroscopic propagator 
l/(— icj + DooQ^) as g 0. The contribution of the even 
terms, n = 2, 4, 6, . . . , is 0{q^) and does not affect the dif- 
fusion coefficient. We now focus on the 0{q'^) contribution 
from the odd n terms: 

1 
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FIG. 6. The asymptotic power-law tail in the time-dependent dif- 
fusion in muscle fibers from ref. 16. (A), Comparing the t~^/^ 
(thin blue and red dashed lines) and (thick black dashed lines) 
power laws for calf tongue genioglossus (TG) and heart (H), cf. 
Fig. 3b. (B), Three-parameter fit of (t) to the power law (5), ac- 
counting for experimental error bars, for tongue (four longest time 
points) and for heart (all data points) yields power law exponents 
'd = 0.44 dz 0.30 and 0.61 ± 0.07 correspondingly. Large standard 
error in tongue is due to a relatively narrow range of experimentally 
available t. 



yielding V{uj) = + cj^F(cj), where 



2C 



1 



= 1,3,. 



1,2 



5(1 + 



kna cot kna 



-luo ■ 



Dokl 



and tan(/i;„a/2) = —Qkna/2 since (/ — > 0. Using equation 
(8), we obtain Anst(t) = -Doo - dtF{t), such that 



2C^ exp(-t/t„) 
•C+(Cfcna/2)2' 



- = Dokl 
(29) 

This expression is exact for all t. For long t>tD= o? /2Dq, 
where is the time to diffuse across one interval, already 
the first term suffices (the others decay exponentially faster), 
agreeing very well with the numerical result in Fig. Ic. 

We note that the residence time in one box, = C'^d, does 
not enter equation (29), so that the times tn are instead deter- 
mined by to- This is remarkable, as it is that is the time 
scale that determines the physics of transport at large tempo- 
ral and spatial scales for any disordered case considered above 
(and Tr ^ to for weakly permeable barriers, C ^ 1). This 
observation suggests that perfectly ordered systems are excep- 
tional, in the sense that the time dependence of transport in 
them is not representative of most biological and man-made 
samples that are at least somewhat disordered. Diffusion (and 
transport in general) in perfectly ordered samples exhibits co- 
herence due to the infinitely long spatial correlations. Physi- 
cally, for the system to equilibrate (establish a long-term den- 
sity profile) it is enough for the density of random walkers to 
equilibrate within each identical interval; there is no need to 
hop over barriers to sense the density in neighboring intervals. 



Permeability and cell size in muscle fibers 

Here we focus on the eigenvalues A2 (t) and Xsit) [notation 
of ref. 16] transverse to muscle fibers in the heart (H) and 
in the tongue genioglossus (TG). (We do not consider here 
the case of tongue core also measured in ref. 16, where the 
geometry is more complicated than that of parallel fibers.) 

As diffusion is axially symmetric [16], X2{t) ^ Xs{t) both 
in H and TG, we consider the isotropic transverse component 
A± = (A2 + A3)/2, contrasted in Fig. 3a with the practically 
nondispersive eigenvalues Ai along the fibers. 

Fig. 6 demonstrates that the asymptotic t ^ 00 behavior 
of X±{t) is consistent with equation (4) with ^ = 1/2. (The 
analysis of A2 (t) and A3 {t) separately, not shown here, yields 
similar results.) In particular. Fig. 6a demonstrates that the 
competing possibilities of 1? > 1 are not consistent with the 
data; ^ = 1 would correspond to A^ {t) Doo + const • 
(In t)/t according to equation (4), whereas ^ > 1 would result 
in X±{t) Doo + const/t according to equation (5). 

In this way, the dynamical exponent (3) helps identify the 
ds = 1 restrictions to the d = 2 dimensional diffusion. To 
quantify the physical parameters of the underlying muscle 
fiber membrane, one needs a model for the time dependent 
diffusion beyond the asymptotic long-time regime. Recently, 
we found the corresponding dispersive diffusivity for the ran- 
dom barriers (hyperplanes) in any dimension d\ 
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(30) 



ref. 15, approximately for all uj, permeability k, and surface- 
to- volume ratio S/V. Its low-frequency behavior is indeed 
characterized by the exponent ^ = 1/2 in any dimension- 
ality d. Here, I^o is the unrestricted diffusion coefficient, 
C = Si/Vd [cf. equation (12)], 2£ = Dq/k, = iViuJr, 
and r = i'^/Do = Do/{2k,)'^. The corresponding expression 
for D{t) is obtained using numerical integration procedure de- 
scribed in ref. 15. 

This allows us to move one step further and quantify the 
membrane permeability and surface-to- volume ratio. The fit 
of the in-plane diffusivity Xj_{t) to the time dependent dif- 
fusion coefficient D{t) obtained from equation (30) with the 
dimensionality d = 2 yields the parameter values summarized 
in the Supplementary Table 1 . 

Tongue genioglossus (TG): the fitted value of the unre- 
stricted diffusivity I^o agrees well with the longitudinal eigen- 
value Ai. This already indicates that the transverse diffusion 
is predominantly restricted by the fiber walls, whereas the 
diffusion within the fibers (at short times) is approximately 
isotropic. 

Heart (H): since the transverse eigenvalues for the heart 
muscle exhibit the dependence (5) for the whole time range, 
for a stable fit we needed to fix one of the fit parameters; we 
chose to set I^o for the heart to the value of the corresponding 
non-time-dependent eigenvalue Ai, assuming that to corre- 
spond to the unrestricted diffusion coefficient within the fibers 
in analogy with the above case of the tongue muscle. 
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TABLE I. Fit results for transverse eigenvalues X±{t) in heart (H) and tongue genioglossus (TG) muscles [16], Fig. 3. First three columns are 
the fit parameters using D{t) obtained from equation (30), the rest are the quantities derived from them. 
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r, ms 


Doo, /xm^/ms 


TD, ms 


Tr, ms 


i, /xm 
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K X 10^, /xm/ms 


a, /xm 


H 


0.8 (fixed) 


2.59 


68.8 


0.223 


20.5 


26.5 


7.42 


0.699 


5.39 


5.72 


TG 


1.09 


2.84 


1328 


0.283 


329 


467 


38.0 


0.150 


1.43 


26.7 



From the surface-to-volume ratio we estimate the typical 
"cell diameter" a 2d/{S/V) which agrees with histological 
values for the actual muscle fiber diameters in both kinds of 
muscle (a = 20 — 40 /xm for tongue and a = 6 — 12 //m for 
heart, see ref. 16 for references), yielding the heart fibers to be 
much narrower. The diameter values are closer to histology 
than those determined in ref. 16 using a fit to an empirical 
two-compartment model. 

Furthermore, our model (30) allows us to determine the 
membrane permeability values k ^ 10~^ cm/s. The k value 
for the tongue muscle agrees well with that expected for cell 
plasma membranes in eukaryotic cells [40]. Unfortunately, 
there is no "gold standard" noninvasive method to determine 
membrane permeability, which makes precise validation of 
the values of n currently unfeasible. However, qualitatively 
and quantitatively, these values are meaningful. In particular, 
the apparent permeability in the heart is almost as large as that 
of a red blood cell membrane [41], a few times more perme- 
able than that of the tongue fibers. Such an elevated value is to 
be expected, since in the heart, there is an abundance of highly 
permeable blood capillaries of a similar diameter aligned with 
fibers, which is likely to increase the average permeability of 
all barriers. We also note that, while our model [15] does not 
include extracellular space, its effect is arguably not crucial, 
as most of the water in tissues (80% or more [42]) is con- 
tained inside cells. Based on the permeability values, on the 
measured time scales water molecules are able to enter and 
exit muscle fibers, rather than being confined within them or 
within the extracellular space, further justifying the use of the 
relatively simple random-membrane geometry of (30). 



Origins of time-dependent diffusion in brain 

Here we demonstrate how the low-frequency dispersion 
(6) with ^ = 1 /2 allows one to radically narrow down the 
scope of the plausible scenarios of the structural or functional 
changes in ischemic stroke. We note that the relatively high 
frequencies (kHz) for dMRI in Fig. 4 appear to be low from 
the point of the brain microstructure at the fim scale, with the 
clinical dMRI accessing D^o = V{uj)\^=o only. 

(i) Active or passive transport? While active axonal trans- 
port, cytoplasmic streaming and microcirculation have been 
discussed as possible reasons for why diffusion might be en- 
hanced in normal tissue relative to postmortem (see recent 
ref. 43 for a review), it seems unlikely that the combination 
of these effects alone could yield the power-law dispersion (6) 
even in a normal case. Indeed, the lack of a time scale in equa- 



tion (6) means that these streaming processes must happen on 
multiple time scales, fine-tuned in such a way as to produce 
the exact power law exponent 1/2. Even if this were the case, 
such fine-tuning must break down after cell death with those 
processes switched off, causing the dispersion to change qual- 
itatively, which contradicts Fig. 4b. Hence we conclude that 
water motion is mostly determined by the ordinary diffusion 
hindered by the passive restrictions, and the change in their 
properties after injury is only quantitative, not qualitative. 

(ii) Order or disorder? The dispersion (6) is non-analytic 
in frequency. Hence, neither bounded motion (water confined 
e.g. to impermeable cells of finite volume, or to effectively 
disconnected pockets of the extracellular space), nor any pe- 
riodic structures (e.g. periodic permeable barriers, periodic 
beads [27], or any packing with a single pronounced length 
scale in any dimensionality [31]) provide the dominant cause 
for the observed dispersion. Indeed, all of those cases cor- 
respond to = oo. The quantity V{uj) in this case is an 
analytic function of co, i.e. it can be Taylor-expanded for 
small uj. As the velocity autocorrelator {v{t)v{0)) is real, 
the measured real part [31, 32] of its Fourier transform must 
be an even function of ou. This means that its Taylor ex- 
pansion at small uj must start with an ou'^ term, RqV{uj) = 
Doo + const • cj^ + 0{uj^). This analytic behavior was ex- 
perimentally demonstrated in ref. 44 for porous samples with 
impermeable walls, and is inconsistent with Fig. 4. Hence, the 
predominant restrictions to diffusion in ref. 20 are nonconfin- 
ing and disordered. 

(iii) Which disorder class? The passive restrictions, while 
not completely ordered, are still correlated in space so as to 
yield the o;^/^ behavior. This observation is crucial: in the 
absence of any structural correlations (p = 0), from equation 
(3) one expects the cj^/^ dispersion in a random d = 3 di- 
mensional medium. Hence, either the disorder is long-range 
correlated (cf. Fig. 2), or the effective dimensionality is less 
than 3. We also note that the nonzero value of rules out 
the "anomalous" diffusion [2], e.g. in a fractal geometry — in 
other words, the structural correlations are gradually forgotten 
and the dynamics asymptotically becomes IMarkoffian. 

The power law exponent ^ = 1/2 can arise due to the two 
remaining classes of the passively restricted diffusion: 

• d > 1: The extended disorder with tig = 1 in = 2 or 
ds = 2 in d = 3, Fig. 2. The permeable barriers may corre- 
spond to either plasma membranes of neurons, glial cells, and 
of their processes (conceptually similar to the above exam- 
ple of muscle fiber membranes), or the membranes surround- 
ing intracelluar organelles, such as nuclear envelope or endo- 
plasmic reticulum. The analysis, based on equation (30), ren- 
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ders this possibility unlikely based both on the corresponding 
length scales and the permeability values. 

• d = 1: Any short-ranged disorder, with the effective di- 
mensionality d = 1, such as the p = example in Fig. 1. 
Below we will argue that this is consistent with an effectively 
one-dimensional water motion along locally straight narrow 
neurites (dendrites and axons) as well as the processes of 
glial cells, assuming their walls to be impermeable, with some 
structural disorder (e.g. beads and shafts) along the way. 

(iv) Extracellular water is less important. The much- 
debated contribution of the extracellular water [45] does not 
contribute to the observed "3=1/2 dispersion. Indeed, its 
effective dimensionality, assuming no exchange with cells, 
would be either d = 2 (due to tight cell packing [42]) or 
d = 3. Any short range disorder (p = 0) in the extracellu- 
lar space would then result in ^ = 1 or 3/2 correspondingly, 
making this contribution less relevant for the observed disper- 
sion: in a superposition of the cj^/^ contribution, and of the 
Ic^l^ or cj^/^ contributions, the cj^/^ dominates as ^ 0. 
This complements the dMRI measurements of various intra- 
cellular metabolites [46-50], indicating that major ischemia- 
related changes occur already in the intracellular space. 

The case of p = and d = 1 (neurites). Let us now 
assume that the neurites are impermeable and locally straight 
narrow one-dimensional channels with water volume fraction 
(j)id, similar to those suggested in refs. 25 and 26. In contrast 
to refs. 25 and 26, we would not assume them to be hollow 
cylinders; rather, we allow structural disorder along the chan- 
nels, leading to the dispersive effective one-dimensional diffu- 
sivity Vid{uj) identical for each channel. Then, the measured 
dispersion of diffusion in a particular direction 

V{cu) ^ ^idVid{oo)/3 + (1 - . (31) 

Here the factor 1/3 assumes approximately isotropic direc- 
tional distribution of neurites in gray matter, and De now is 
the effective diffusivity of water outside the neurites. We can 
approximately set De = const, as the residual dispersion in 
the extra-neurite space should be less singular than cj^/^ since 
its dimensionality d > 1 (cf. our discussion in (iv) above). We 
then focus on the one-dimensional channel dispersion Vid{(ju) 
from equation (31), which imposes constraints on the possible 
parameter values. 

Our main constraint will be on the neurite volume fraction 
(pid. First, since ReVid{(ju) > 0, 

(1 - < V^in = V{0) , (32) 

where the long-time limit X>min ~ 0.74//m^/ms in normal 
and Pmin ~ 0.5 fim^/ms in globally ischemic brain. Fig. 4. 
Second, the d = 1 diffusivity ReVid{uj) < D^yt cannot ex- 
ceed that of water in cytoplasm, D^yt < 3/im^/ms; from 
ref. 51 the axoplasm diffusivity is about 80% of that of pure 
water, yielding D^yt ~ 2.4 /im^/ms as a plausible estimate. 
Thus from equations (31) and (32), we obtain 

3(I>max - ^min)/01d < Dcyt , (33) 



where Pmax is the maximal measured ReP(cj), Fig. 4. 
Taking into account the measurements up to IkHz of the 
cos waveform [20], the range is at least X>max — ^min ~ 
0.3 /im^/ms for both normal and globally ischemic brain, 
while the diffusion length at 1 kHz still exceeds the neurite 
inner diameter of a fraction of /im, so that the motion within 
the channel remains one-dimensional. Hence, from equation 
(33) it follows that in order to achieve the observed dispersion 
with one-dimensional neurites, their volume fraction should 
be sufficiently large, (j)id ^ 0.4. This is consistent with the 
neurite volume fraction (j)id ~ 0.6 measured with electron 
microscopy [52]. 

Which physical parameter changes most with ischemia? 
First we note that the observed 50% increase of the coefficient 
in front of cj^/^ under ischemia. Fig. 4, generally signifies a 
relative increase of the disorder, {{SD)'^) / D^^, causing the 
cj^/^ dispersion [cf. (7)]. One consequence of this observa- 
tion is that this apparent disorder increase is inconsistent with 
another debated scenario of why drops in stroke — a sug- 
gestion that the cytoplasm itself becomes more "viscous", or 
"dense", causing the decrease of the free diffusion Do in is- 
chemia, and with that, of I^oo • Physically, the decrease of I^o 
would either not affect the structural disorder or cause its rel- 
ative decrease, by reducing the contrast between regions with 
freely diffusing water and any bottlenecks or barriers. Indeed, 
as a concrete example, the pref actor of cj^/^ in the small-cj ex- 
pansion of equation (30) decreases with the decrease of Do for 
any dimensionality d. Hence, the signature decrease in in 
ischemia cannot be explained by assuming that the cell cyto- 
plasm becomes "denser" or "more viscous" after injury. This 
again hints at the major changes being structural rather than 
molecular. 

A crude estimate for the change in the one-dimensional dis- 
order after ischemia could be made using our model [15] of 
randomly placed barriers, equation (30), ind = 1 dimension. 
For that, we assume that, on the length scales longer than a few 
/im, the narrow shafts between beads and spines act as effec- 
tive barriers. Their effective one-dimensional "permeability" 
would be smaller for either narrower shafts, or thicker beads, 
and grows with the ratio between the diameters of the shafts 
and the beads. 

The difficulty here is that the empirical dependence in Fig. 4 
allows one to obtain only two parameters, as the measurement 
[20] does not extend towards high enough uj for which V{uj) 
saturates. Hence, as above, we need to ensure that the pa- 
rameters such as (l)id and De stay within the above bounds. 
With that in mind, we fit the real part of equation (30) with 
d = 1 to Vid{io) from (31). Choosing (t)id and De within 
the above bounds (so that, e.g., ReVid{u)) < Dcyt), we ob- 
tain reasonable fit results by setting (j)id = 0.7 for both before 
and after global ischemia, this fraction being slightly greater 
than that observed in the neuropil [52]. This requires the dif- 
fusivity De outside the neurites to decrease from about 1.9 
before to 1.2 /im^/ms after ischemia. This may be explained 
by the effect of cell swelling making the geometry more "tor- 
tuous" [42] which, as argued in (iv) above, does not affect 
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the cj^/^ dispersion. We find that the effective barrier perme- 
ability K drops from about 0.5 /xm/ms before to 0.2 /xm/ms af- 
ter ischemia, while the distance between "barriers" increases 
from a = 2 to 3 /im. The permeability reduction is consis- 



tent with the more pronounced beads (stronger contrast be- 
tween bead and shaft diameters) in injured dendrites, and the 
increase in the distance is consistent with the disappearance of 
spines [23], such that the disorder correlation length increases. 



